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Abstract 

In the Minimal Super symmetric Standard Model (MSSM) we incorporate the Higgs- 
boson propagator corrections, evaluated up to two-loop order, into the prediction of 
T(h — > ff) and BR(h — > //) for f = b,c, r. The propagator corrections consist of the 
full one-loop contribution, including the effects of non-vanishing external momentum, 
and corrections of 0(aot s ) at the two-loop level. The results are supplemented with 
the dominant one-loop QED corrections and final state QCD corrections from both 
gluons and gluinos. The effects of the two-loop propagator corrections and of the one- 
loop gluino contributions are investigated in detail. Our results are compared with the 
result obtained within the renormalization group approach. Agreement within 10% is 
found for most parts of the MSSM parameter space. 
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1 Introduction 



The search for the lightest Higgs boson is a crucial test of Supersymmetry (SUSY) that can 
be performed with the present and the next generation of accelerators. The prediction of a 
relatively light Higgs boson is common to all supersymmetric models whose couplings remain 
in the perturbative regime up to a very high energy scale |]J . Finding the Higgs boson is thus 
one of the main goals of todays high-energy physics. Concerning the Higgs boson search, it is 
necessary to know the decay widths and branching ratios of the main decay channels to a high 
accuracy. After the detection of a scalar particle it is mandatory as a next step to measure 
its couplings to gauge bosons and fermions and also its self-couplings very accurately, in 
order to establish the Higgs mechanism and the Yukawa interactions experimentally. The 
determination of the trilinear Higgs-boson self-couplings might be possible at a future linear 
e + e~ collider with high luminosity M. 

In this paper we concentrate on the coupling of the lightest MSSM Higgs boson to 
Standard model (SM) fermions. In the MSSM the mass of the lightest Higgs boson, m^, 
is bounded from above by rrih ^$ 135 GeV, including radiative corrections up to two-loop 
order [[3-11]. Since the b-, the c-quark and the r-lepton are the heaviest particles for which 
the decay h — > ff is kinematically allowed, it is of particular interest to calculate their 
corresponding decay rates and branching ratios with high precision 

EMI- We analyze 



these decay rates and branching ratios, taking into account the Higgs-boson propagator 
corrections, where at the one- loop level the full momentum dependence is kept. These 
corrections contain the Yukawa contributions of Oi y Gpm\jM'^ v \ which are the dominant 
electroweak one-loop corrections to the Higgs-boson decay width, and the corresponding 
QCD corrections of 0{ y Gpa s m\lMy V ) as well as the Yukawa corrections of 0{G 2 F m^ /M^). 
We also take into account the one-loop vertex corrections resulting from gluon, gluino and 
photon exchange together with real gluon and photon emission as given in Ref. |12| . Only 
the purely weak 0(a) vertex corrections have been neglected. We numerically investigate 
the effect of the two-loop propagator contributions and the one-loop gluino-exchange vertex 
correction. The latter one has been mostly neglected in experimental analyses so far, but can 
have a large impact on the result. We show analytically that the Higgs-boson propagator 
correction with neglected momentum dependence can be absorbed into the tree-level coupling 
using the effective mixing angle from the neutral CP-even Higgs boson sector. The result 
in this approximation is then compared with the full result. We also perform a comparison 
between our full result and the renormalization group (RG) improved effective potential 
calculation, see Ref. [T9I] and references therein. For most parts of the MSSM parameter 



space we find agreement within 10% between the two approaches, although deviations up to 
50% are possible for certain ranges of the parameter space. 

The paper is organized as follows: in section 2 we give the calculational basis needed 
for the incorporation of the two-loop Higgs-propagator corrections into the decay widths. 
We show analytically how the propagator corrections with neglected external momentum 
are related to the effective mixing-angle approach. The results for the gluon, gluino and 
QED vertex corrections in combination with gluon and photon bremsstrahlung are reviewed. 
In section 3 a numerical analysis for the decay rates and a comparison of the full result 
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and the effective mixing-angle result is performed. Special emphasis is put on the two- 
loop propagator correction and the one-loop gluino contribution. Section 4 contains the 
comparison of our full result with the RG approach. The conclusions can be found in 
section 5. 



2 Calculational basis 

The Higgs potential of the MSSM is given by 

V = m\H 1 H 1 + mlH 2 H 2 -m\ 2 {e ab H^H b 2 + h.c.) 

+ 9 —^- {H l H l - E 2 E 2 f + £- IH.H.I 2 , (1) 
o Z 

where m 1 ,m 2 ,m 12 are soft SUSY-breaking terms, g,g' are the 577(2) and £7(1) gauge cou- 
plings, and ei2 = —1. The doublet fields H 1 and H 2 are decomposed in the following way: 

H ( H l\ - ( vi + W + il&yV* \ 

H2 = (^) = (^ + (^+<xS)m)- (2) 

Besides g and g', two independent parameters enter the potential fll]): tan/3 = v 2 /vi and 
M\ = — (tan /3 + cot (3 ), where Ma is the mass of the CP-odd A boson. 

The CP-even neutral mass eigenstates are obtained performing the rotation 
/ H° \ ( cos a sin a \ ( (p\ \ _ n _i, >. / 4>i 



h° 



o 



— sin a cos a I \ <P 2 I \ <P2 



D ~w 33 ( 3 ) 



with the mixing angle a related to tan (3 and Ma by 

Ml + Ml 7T 

tan2a = tan2/3 M ^_ M f , -- < a < 0. (4) 

At the tree level the mass matrix of the neutral CP-even Higgs bosons in the <fri,<f) 2 basis 
can be expressed in terms of Mz and Ma as follows: 



Mi 



2,tree _ / m\ 



Higgs I 2 2 



Ml sin 2 /3 + Mf cos 2 /3 — (M| + Mf) sin /3 cos (3 
-(Ml + Ml) sin (3 cos /3 Mj cos 2 (3 + M| sin 2 /3 



(5) 



Transforming to the eigenstate basis (^) yields 



Mgr^M m n * \ i ((») 



with rrih and m# being the tree-level masses of the neutral CP-even Higgs bosons. 



In the Feynman diagrammatic (FD) approach the higher-order corrected Higgs boson 
masses, denoted by M h , M H , are derived by finding the poles of the h, if -propagator matrix 
whose inverse is given by 

(A r l __■( q 2 -m 2 H + t H (q 2 ) t hH (q 2 ) \ 

1 HiggsJ " V W) q 2 -< + Uq 2 ))' {) 

where the S(g 2 ) denote the renormalized Higgs boson self-energies. For these self-energies 
we take the result up to two-loop order, see Sect. below. 

Our main emphasis in this paper is on the fermionic decays of the light Higgs boson, but 
for completeness we list the expressions for both h and H . The amplitudes for the decays 
h,H—>ff can be written as follows: 



with 



A(h - ff) 
A(H -> ff) 

ZhH = 
Zuh = 



Z h (T h + 

ZhH Tj?) 
Zh (Th + Zffh Th) 

t hH (Ml) 



M 2 h -m 2 H + Z H (M* h ) 



t hH (M 2 H ) 



M 2 H -m 2 h + Y, h (M H ) 



(9) 

(10) 
(11) 



involving the renormalized self-energies tl(q 2 ) and 3-point vertex functions Th, Th- The wave 
function renormalization factors Zh and Zh are related to the finite residues of the h and H 
propagators, respectively: 
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l + ReS' h (g 2 ) - Re 
1 



q 2 -m 2 H +H H {q 2 ) 



'H 



l + Ret' H (q 2 ) - Re 



q 2 -ml+T, h (q 2 ) 



(12) 
(13) 



2.1 The a e ff-approximation 

The dominant contributions for the Higgs boson self-energies can be obtained by setting 
q 2 = 0. Approximating the renormalized Higgs boson self-energies by 



t( q 2 ) - E(0) = t 



(14) 



yields the Higgs boson masses by re-diagonalizing the dressed mass matrix 

M 2 (m 2 H -±H -±hH \Aa(M 2 H \ 

Hisgs "l -±hh ml-t h )^\ M 2 J' ^ 

where and Mu are the corresponding higher-order-corrected Higgs boson masses. The 
rotation matrix in the transformation (1151) reads: 



D(Aa) = C ° S ^ ) . (16) 

V smAa cosAa / 

The angle Aa is related to the renormalized self-energies and masses through the eigenvector 
equation 



m\-Y, H -M 2 h -E hH \f -sinAa 
-t hH m 2 h -t h -M?)[ cosAa 



(17) 



which yields 



Ml -m 2 H + H H 
The second eigenvector equation leads to: 



M 2 H -m 2 h + E h 
With the approximation ( |i"4|) one deduces 



tanAa . (18) 



tanAa . (19) 



Z hH = - Sfe f =- tan Aa, (20) 

Z H h = ^ ~ = +tanAa, (21) 

M 2 H -m 2 h + H h 



and Zh can be expressed as 

Z h 



1 + 



1 



cos 2 Aa . (22) 



Analogously one obtains 



1 + tan 2 Aa 

Z H = cos 2 Aa . (23) 
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At the tree level, the vertex functions can be written as 

> for d-type fermions 



rv 



lenit sin a 



2swM\y cos j3 
—ierrif cos a 
2swMyy cos /3 

—ierrij- cos a 
2s\yMw sin/3 
—ierrif sin a 



Cy sina 

- —C^p cos a 

C j COS (X , 

- C f sina 



> for u-type fermions . 



(24) 
(25) 



2s\yMw sin /3 

Incorporating them into the decay amplitude yields: 
A eS (h^ff) = 



A cS (h 
A cS (H 
A eS (H 

Recalling the relations 
and 



Zh (Xh + Z hH T H ) 

= Cf 1 cos Aa (sin a — tan Aa (— cos a )) 

= Cf ] sin(a + Aa) 

= Cp* sin a e g (for d-type fermions) 
= cos a cS (for u-type fermions) 

= — Cj^coscteff (for d-type fermions) 

= Cj"' 1 sin a cff (for u-type fermions) 



D(a eff ) = D(a) D(Aa) 



a 



D(a) 



it is obvious that a e ff = (a + Aa) is exactly the angle that diagonalizes the 
corrected Higgs boson mass matrix in the 0i, 2 -basis: 



(26) 
(27) 
(28) 
(29) 

(30) 
(31) 

higher order 



m 



a 



2 v 



m <Pl<f>2 ^4>l<t>2 



a c ff 



Ml 





Ml 



(32) 



m 2 H -H 



H 



— E 



hH 



Aa 



Ml 



mt 



Ml 



a c s can be obtained from 



a c g = arctan 



-(M\ + Mf) sin/3 cos/3 - £^ 2 



Mf cos 2 /3 + Mi sin 2 - E 01 - m 2 , 



7T 

- - < «cfT < 



7T 

2 



(33) 
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2.2 The Higgs-boson propagator corrections 

For the Higgs boson self-energies employed in eqs. ©-(0) we use the currently most 
accurate result based on Feynman-diagrammatic calculations. It contains the result of 
the complete one- loop on-shell calculation of Ref. [|j], together with the dominant two- 
loop corrections of 0(aa s ) obtained in Refs. fllQjllll, including also the leading terms of 
0{G 2 F m^/ 'M^y) the Fortran program FeynHiggs, based on this result, has been de- 

scribed in Ref. In this way the complete MSSM one-loop on-shell result together with 
the dominant two-loop contribution, originating from the t — t-sector (without any restric- 
tions on the mixing), is taken into account in the propagator corrections. 
In the approach in Refs. [|K,TT] the Higgs boson self-energies are given by: 



± s (q 2 ) = ±^(q 2 ) + ±f\0) , s = h,H,hH, (34) 

where the momentum dependence has been neglected only at the two-loop level, while the 
full momentum dependence is kept in the one-loop contributions. 

In a first step of approximation for the calculation of the decay width T(h — > ff) the mo- 
mentum dependence is neglected everywhere in the Higgs boson self-energies (see eq. (14)): 

± s (q 2 ) -> tW(0) + Si 2) (0) , s = h,H,hH. (35) 



This corresponds to the a e ff-approximation, as described in Sect. 2.1 



In a second step of approximation we approximate the Higgs boson self-energies by the 



compact analytical formulas given in Ref. |p2 



± s (q 2 ) = S(!)™ x (0) + S(2)approx( ) ? S = h,H, hH, (36) 

yielding relatively short expressions which allow a very fast numerical evaluation. In the 
following, this approximation is labeled by a e ff(approx). 



For the t-sector, we use the same conventions as in Ref. 11]: the scalar top masses and 



the mixing angle are related to the parameters M$ , Afj and X t of the t-mass matrix 

KA* - ( M t + m t + C ° S2 ^ (5 - ¥w) M2 Z m ^ \ (a* 

JV[ i-{ m t X t M? R +m 2 t + 1 cos2(3 s 2 w M 2 ) ' 

with 

X t = A t - fx cot p . (38) 

In the numerical analysis below we have chosen = Mi L = M^ R . 

In Refs. it has been shown that for a given set of MSSM parameters the maximal 

values of Mh as a function of X t are obtained for \X t /mq\ w 2. This case we refer to as 
'maximal mixing'. Minimal values for Mh are reached for X t ~ 0. This case we refer to as 
'no mixing'. 
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2.3 Decay width of the lightest Higgs boson 



At the tree level, the decay width for h — > ff is given by 



r (^ - //) = n< 



m h 



1 - 



4m 2 , 



. 7T 



mi 



(39) 



The electroweak propagator corrections are incorporated by using the higher-order decay 
amplitude (0) 



The a e fr-approximation is given by 



iV, 



1 - 



4m 



. 7T 



M, 2 



l^-//)| 5 



eff 



> 7T 



(40) 



(41) 



In this paper we consider only those electroweak higher-order contributions which enter 
via the Higgs boson self-energies. These corrections contain the Yukawa contributions of 
0(GFm'f/M^ v ), which are the dominant electroweak one-loop corrections to the Higgs-boson 



decay width, and the corresponding dominant two-loop corrections, see Sect. 2.2. The pure 



weak 0(a) vertex corrections are neglected (they have been calculated in Ref. [0 and were 
found to be at the level of only a few % for most parts of the MSSM parameter space, see 



also Sect. |2.3.3| below.) 



2.3.1 QED corrections 



Here we follow the results given in Refs. [12-15]. The IR-divergent virtual photon contri- 
bution is taken into account in combination with real-photon bremsstrahlung yielding the 
QED corrections. The contribution to the decay width induced by 7-exchange and final-state 
photon radiation can be cast into the very compact formula 



Ar 7 = r\ • 5T y , 
where for m 2 <C M| the factor 5Tj has the simple form 

a 



6T„ 



3 log m + 1 

\m f J 4 



(42) 



(43) 



2.3.2 QCD corrections: gluon contributions 

The corresponding results have been obtained in Refs. flT2"|-|nfj. The additional contribu- 
tion to the decay width induced by gluon exchange and final-state gluon radiation can be 
incorporated into (0) by writing 



m 2 JMl 



mi 




1 - 



(44) 



7 



The correction factor containing a s (m 2 ), which has not been included in previous diagram- 
matic calculations, can give rise to non-negligible contributions. m q (M%) is calculated via 



m q (q 2 



m 



c{q 2 ) 
q c(m q ) 

' P a s (q 2 
2tt 



-7o/2/3 



1 + 



(Pijo ~ Poll) a s {q 2 



Pi 



8tt 



+ 



' {Pilo - Polif 7o(p2po - Pi) JlPl _ J2 

Po Po Po , 



2/3 4 



> 7T 



(45) 



(46) 



where m q is the pole mass and m q {m q ) = m q . The coefficients in eq. 

33 - 2N f 



Po 

Pi 

Pi 
To 

7i 

72 



102 - —Nf 

3 ; 

2857 5033 



2 

•8 , 
404 



-N f + 



325 
5? 



-N 



f > 



40 *r 

1 Nf , 

3 9 ' ' 

^Nj + (l60C(3) + ™)N f - 3747 



61) are: 



(47) 



where ((3) ~ 1.2020596 . . . and Nf = 5 for f = c,b, which are considered here. The strong 
coupling constant a s is given up to three loops by: 



An 
PoL q 



_ P^logLq P 2 log 2 L q _ /glogL, p 2 p - P 2 1 

- * tr%A t O /nA t- O ' 



Pt Ll 



P'o LI 



L\ 



(4S 



where L q = \og(q 2 / Aq CD ) . (For the numerical evaluation Aq CD = 220 MeV has been used.) 
Numerically, more than 80% of the gluon-exchange contribution is absorbed into the running 
quark mass. 



2.3.3 QCD corrections: gluino contributions 

We follow the calculation given in Ref. \\L2jf\ similar results can also be found in Ref. fl~7j| . 
The additional contributions to the decay width induced by gluino-exchange are incorporated 
via 



AT; 



Ti ■ 5T- g , 



where STg is given by 



6T S 



T/i + Z hH Y H 



-Re 



It + Z h 



(49) 
(50) 



1 An error in Ref. ]12|] concerning the proper inclusion of the Hff coupling has been corrected. 
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for real Z^h (i.e. neglecting the imaginary part in eq. fl50|)); r~ and are given by 



AT~ h 



1 2 =M 2 h 



+ Ei ^ ( m 



/' / 2 



r 



2- 



(51) 



(52) 



AT^ ,H denote the gluino vertex-corrections, whereas represents the gluino contribution 
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to the fermion self-energy corrections. Explicit expressions for these terms can be found 



in Ref. 12 



For large values of tan/3 in combination with large values of the gluino-exchange 
corrections to T(h — > bb) can become very large. In Refs. [[23|,|24| as well as in Ref. [25| it 
has been proposed to derive an effective contribution to the decay width resummed to all 
orders. A similar resummation can be applied for weak 0(a) chargino-exchange corrections 
to Y(h — > bt>), where they become non-negligible ||23|| . In the numerical examples given 
in Ref. |25| the difference between the one-loop result for T(h — > //) and the effectively 
resummed result does not exceed 10 - 15%, even for very large values of tan/3 > 40. A 
proof of how the resummation of the leading terms arising for large tan f3 can be performed 
for the H + ib vertex is given in Ref. [26|. Concerning our numerical analysis in Sect. |] 



we have neglected the additional contributions from the resummation. These additional 
corrections, although potentally large in some regions of the MSSM parameter space, would 
not qualitatively change our conclusions given below. A more detailed investigation of the 
effects of a proper resummation of the leading contributions and of the inclusion of the 
complete electroweak one-loop vertex corrections will be given in a forthcoming publication. 
For some parameter combinations the gluino corrections can drive T(h — > bb) to very 



small values, see the discussion at the end of Sect. 3.1 



2.3.4 Decay width and branching ratio 

Including the various types of corrections, the decay width is given by 

r(/i - //) = r li9 + Ar 7 + at § . (53) 



Summing over / = 6, c, r and adding T(h — > gg) (which can be numerically relevant |[27||), 
results in an approximation for the total decay width 

r to t= E nh^ff)+T(h^gg) . (54) 

f=b,c,r 

We do not take into account the decay h — > AA (see e.g. Ref. for a detailed study). 
Although it is dominant whenever it is kinematically allowed, it plays a role only for very 
small values of tan/? (tan/? < 1.5) which will not be considered here because of the limits 
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obtained at LEP2 [29]. We also assume that all other SUSY particles are too heavy to allow 
further decay channels. In addition, we neglect the decay h — > WW* which can be of 0(1%) 
for M h > 100 GeV. 

The fermionic branching ratio is defined by 

R, = BR{h -> //) = F(/l T T /7) • (55) 

1 tot 

3 Numerical analysis 

Concerning the numerical evaluation of the Higgs-boson propagator corrections, we follow 



Sect. |27|. For tan (3 we have chosen two representative values, a relatively low value, tan {3 = 

3.0 and a high value, tan/? = 40. For sake of comparison we also consider an intermediate 
value of tan f3 = 20 in some cases. If not indicated differently, the other MSSM parameters 
are chosen as follows: \i = —100 GeV, M 2 = rriq (M 2 is the soft SUSY-breaking term 
in the gaugino sector), gluino mass m~ g = 500 GeV, Aj, = A t (which fixes, together with 
fi, the mixing in the 6-sector). For the SM fermion masses we have furthermore chosen 
m t = 175 GeV, m b = 4.5 GeVfl , m T = 1.777 GeV and m c = 1.5 GeV. 

The mass Ma of the CP-odd Higgs boson is treated as an input parameter and is varied 
in the interval 50 GeV < Ma < 500 GeV. The corresponding values for Mh follow from 
eq. (^). Mh, derived in this way, subsequently enters the numerical evaluation of the formulas 
presented in section Thus the variation of Mh in the plots stems from the variation of Ma 
in the above given range. 

3.1 Effects of the two- loop Higgs-propagator corrections 

We first focus on the effects of the two-loop Higgs-boson propagator corrections. They have 
been evaluated at the one- and at the two- loop level as described in Sect. |27|. Fig. [I] shows 
the results for T(h — > bb) for a common scalar quark mass rriq = 1000 GeV and tan/3 = 3 
and tan /3 = 40 in the no-mixing and the maximal-mixing scenario. The QED and the QCD 
gluon and gluino vertex contributions are also included. 

In the small tan/5 scenario, larger values for T(h — > bb) are obtained for maximal mixing. 
The two-loop corrections strongly reduce the decay width. In the large tan/5 scenario the 
variation is mainly a kinematical effect from the different values of Mh at the one- and two- 
loop level. The absolute values obtained for T{h — > bb) are three orders of magnitude higher 
in the tan/5 = 40 scenario, which is due to the fact that T(h — > bb) ~ 1/ cos 2 /? . 

In Fig. [| the three decay rates T{h — > bb), Y(h — > t + t~) and T(h — > cc) are shown as 
a function of Mh- The results are given in the no- mixing scenario for rriq = 500 GeV and 

2 This values is well above the expected limit obtainable at the end of LEP, assuming that no Higgs boson 
signal will be found p3]. For these expected limits m t — 174.3 GeV and rriq — 1 TeV has been assumed. 
The value of rrih = 4.5 GeV is used in the tree-level expression and in the QED and QCD vertex 
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corrections (see Sects. 2.3.1 - 2.3.3), while for the Higgs-propagator corrections the running bottom mass, 



Ti(,(mt) = 2.97 GeV, has been used, in order to partially absorb higher-order QCD corrections. 
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Figure 1: T(h — > bb) is shown as a function of Mh- The Higgs-propagator corrections have 
been evaluated at the one- and at the two-loop level. The QED, gluon and gluino contributions 
are included. The other parameters are ji = —100 GeV, M 2 = mq, rrig = 500 GeV, A b = A t , 
tan p = 3,40. The result is given in the no-mixing and maximal-mixing scenario. 

tan/5 = 3,40. 

In the low tan/? scenario T(h — > bb) and T(h — > t + t~) are lowered at the two-loop level, 
while T(h — > cc) is increased. The decay rate for h — > bb is about one and two orders of 
magnitude larger compared to the ones of h — > t + t~ and h — > cc, respectively. In the large 
tan/3 scenario the shifts are again dominated by the kinematical effect from the different 
values of Mh at the one- and two-loop level. In the maximal-mixing case, which is not 
plotted here, we find qualitatively the same behavior. 

We now turn to the effects of the two-loop corrections to the branching ratios. For not too 
large values of Mh, T tot is strongly dominated by T(h — > bb). For large values of Mh the decay 
into gluons becomes more relevant. In Fig. [| we show the branching ratio BR(h — > bb) as 
a function of Mh and Ma- For values of Ma ^ 250 GeV there is a non- negligible difference 
between one-loop and two-loop order, where at the two-loop level the branching ratio is 
slightly enhanced. Compared in terms of Mh there is nearly no change for small values of 
Mh in the low and in the high tan f3 case. Here BR(h — > bb) is changed by less than about 
1%, see Fig. [3|. BR{h — > r + r _ ) is increased by less than about 2%. BR{h — > cc) can be 
increased at the two-loop level by 0(50%), but remains numerically relatively small. For 
tan /3 — 40 the main difference arises at the endpoints of the spectrum, again due to the fact 
that different Higgs boson masses can be obtained at the one-loop and at the two-loop level. 
For tan /3 — 3, however, also several GeV below the kinematical endpoints there is a sizable 
effect on BR(h — > bb). Thus, in the experimentally allowed region of Mh, the two-loop 
corrections can have an important effect on BR(h — > bb). 
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Figure 2: T(h — > bb), T(h — > t + t~) and T(h — > cc) are shown as a function of M^. The 
Higgs-propagator corrections have been evaluated at the one- and at the two-loop level. The 
QED, gluon and gluino contributions are included. The other parameters are /i = —100 GeV, 
M.2 = trig, trig = 500 GeV, Ab = A t , tan/3 = 3,40. The result is given in the no-mixing 
scenario. 

Higgs boson search, especially at e + e~ colliders, often relies on b search, since on one hand 
the lightest CP-even Higgs boson decays dominantly into bb and on the other hand b tagging 
can be performed with high efficiency For some combinations of parameters, however, 
r(/i — > bb) can become very small and thus BR(h — *> bb) can approach zero as a consequence 
of large Higgs-boson propagator corrections or large gluino vertex-corrections, making Higgs 
boson search possibly very difficult for these parameters. Within the effective potential 
approach this kind of effect has first been observed in Ref . p0[ , recent analyses investigating 



the parameter regions where BR(h — > bb) is suppressed can be found in Refs. [23,24]. In 
order to have reliable predictions for these regions of parameter space a full calculation 
of the one-loop vertex corrections, including all 0(a) contributions, would be necessary. 
Here we demonstrate the effect of the two-loop propagator corrections on the values of the 
parameters, especially of Ma, for which BR(h — > bb) goes to zero. We also show the impact 
of the inclusion of the momentum dependence of the Higgs boson self-energies (see eq. (0) 
and (|j~4"D), that is often neglected in phenomenological analyses of the decays of the lightest 
CP-even Higgs boson. 

In Fig. |] BR(h — > bb) is shown as a function of Ma- The Higgs boson self-energies are 
evaluated at the one-loop and at the two-loop level with and without momentum dependence 
(see eq. (|I4"D). The other parameters are tan/3 = 25, rriq = 500 GeV, m g = 400 GeV, 
M 2 = 400 GeV, X t = 400 GeV, A b = A t , fi = -1000 GeV. The inclusion of the two- 
loop propagator corrections shifts the Ma value for which BR(h — > bb) becomes very small 
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Figure 3: BR{h — > bb) is shown as a function of Ma and M h for the same settings as in 
Fig. The QED, gluon and gluino contributions are included. 



by about —35 GeV. The inclusion of the momentum dependence of the Higgs boson self- 
energies induces another shift of about —6 GeV. In order to have reliable phenomenological 
predictions for the problematic Ma values the two-loop corrections as well as the inclusion 
of the momentum dependence is necessary. Note that the inclusion of the gluino vertex 
corrections as well as the purely weak vertex corrections can also have a large impact on the 
critical Ma values. 
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Figure 4: BR(h — > bb) is shown as a function of The Higgs boson self-energies are 

evaluated at the one-loop and at the two-loop level with and without momentum dependence 
(see eq. The other parameters are tan/3 = 25, nig = 500 GeV, m g = 400 GeV, 

M 2 = 400 GeV, X t = 400 GeV, A b = A t ,u, = -1000 GeV. 

3.2 Effects of the gluino vertex corrections 

In this subsection we present the effect of the gluino-exchange contribution to the hff 
vertex corrections. These corrections have been neglected so far in most phenomenological 
analyses^. 

Fig. |5| shows T(h — > bb) in three steps of accuracy: the dotted curves contain only 
the pure self-energy correction, the dashed curves contain in addition the QED and the 
gluon-exchange correction. The solid curves show the full results, including also the gluino- 
exchange correction. The results are shown for the no mixing scenario, \i = —100 GeV, 
nig = 500 GeV, nig = 200, 1000 GeV in the left and right part of Fig. [5], respectively. For 
tan ft, three values have been chosen: tan ft = 3, 20, 40. 

The left plot of Fig. [5| corresponds to a small soft SUSY-breaking scale, nig = 200 GeV. 
The effect of the gluon contribution is large and negative, the effect of the QED correction 
is small. For this combination of m g , nig and \i the effect of the gluino correction is large 
and positive as can be seen from the transition from the dashed to the solid curves. For 
tan/3 = 40 it nearly compensates the gluon effect, for tan/3 = 20 it amounts up to 20% 
of the gluonic correction, while for tan ft — 3 the gluino-exchange contribution is negligible. 

4 The gluino-exchange contributions are currently incorporated into HDECAY Concerning the 

discovery potential of LEP, the Tevatron and the LHC, the gluino corrections have also been studied recently 
in Refs. 0S. 
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Figure 5: F(h — > bb) is shown as a function of M h for three values of tan (3. The Higgs- 
propagator corrections have been evaluated at the two-loop level in the no-mixing scenario. 
The dotted curves shows the results containing only the pure self- energy corrections. The 
results given in the dashed curves in addition contain the QED correction and the gluon- 
exchange contribution. The solid curves show the full result, including also the gluino cor- 
rection. The other parameters are fj, = —100 GeV, M 2 = rriq, Aj, — A t . 
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Figure 6: F(h — > bb) is shown as a function of ji for tan/3 = 20 and two different values 
of trig in the no-mixing and the maximal-mixing scenario. The other parameters are rriq = 
500 GeV, M 2 = 500 GeV, A b = A t , M A = 100 GeV. 
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Note that we have chosen a relatively small value of fi, fj, = —100 GeV. For larger values of 
\fj,\ even larger correction can be obtained. Hence neglecting the gluino-exchange correction 
in the large tan/3 scenario can lead to results which deviate by 50% from the full 0(a s ) 
calculation (see also Sect. |2.3.3| ). The right plot of Fig. |5] corresponds to rriq = 1000 GeV. 



The gluino-exchange effects are still visible, but much smaller than for ms = 200 GeV. The 



same observation has already been made in Ref. |12| . In the maximal- mixing scenario we 
find qualitatively the same behavior for the gluino-exchange corrections as in the no-mixing 
scenario. 

In Fig. |B| the pure gluino-exchange effect is shown as a function of fi. This effect increases 
with risingQmg and |//|, where for negative (positive) /i there is an enhancement (a decrease) 
in T(h — > bb). The size of the gluino-exchange contribution also depends on Ma, where 
larger effects correspond to smaller values of Ma, see also Ref. fll^f . The small difference 
between the curves where the decay rate has been calculated without gluino contribution is 
due to the variation of Mh induced by different values of m g which enters at (D(aa s ). Fig. |] 
demonstrates again that neglecting the gluino contribution in the fermion decay rates can 
yield (strongly) misleading results. 

The gluino exchange contribution has only a relatively small impact on BR(h — > bb). It 
can have a large influence, on the other hand, on BR(h — > t + t~). Both branching ratios are 
expected to be measurable at the same level of accuracy, see e.g. Ref. [ft^] . While the Higgs- 
propagator contributions are universal corrections that affect T(h — > bb) and T(h — > t + t~) 
in the same way (i.e. the influence on the effective coupling is the same in both cases), the 
gluino corrections, which influence only T(h — > bb), can lead to a different behavior of the 
two decay widths. In Fig. |^ we show BR(h — > r + r~) as a function of m g . The left plot 
corresponds to three different values of tan/3 and fi = ±100 GeV, Ma = 100 GeV. Here we 
have furthermore chosen rriq = 500 GeV and moderate mixing, i.e. X t = rriq] we find similar 
results for the other mixing scenarios. For positive (negative) fi the decay rate T(h —>■ bb) is 
reduced (enhanced), see Fig. §, thus BR{h — > t + t~) is increased (decreased) by up to 50%. 
For large values of m g and tan /3, BRih — > t + t~) can thus be considerably different from the 
case where the gluino-exchange contribution to T(h — > bb) has been neglected. This becomes 
even more apparent in the right plot of Fig. 0, where the MSSM result, including gluino 
exchange contribution, is compared to the SM result. Here we also show the scenario with 
Ma = 300 GeV and /i = -400 GeV, where the Higgs sector of the MSSM behaves SM like 
(i.e. the lightest Higgs boson has almost SM couplings, all other Higgs bosons are heavy). 
The horizontal lines represent the SM values for the respective Higgs boson masses. The two 
Higgs boson masses for each line give similar results so that the lines are indistinguishable 
in the plot. These masses correspond to an averaged value obtained in the MSSM in the 
interval < m g < 1000 GeV, where the variation of Mh is about ±1 GeV (with our choice 
of rriq and X t ). Since the gluino decouples very slowly, there is no decoupling effect for 
m g < 1000 GeV and the MSSM results can be considerably different from the SM result, 
for Ma = 100 GeV as well as for M A = 300 GeV, i.e. even where the MSSM Higgs sector 



This is correct for all values of mj considered in this work. A maximal effect is reached around 
» 1500 GeV. The decoupling of the gluino takes place only for very large values, rrig > 5000 GeV. 
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behaves otherwise SM like. The deviation can amount up to 50% for M A 
up to 30% for M A = 0(300 GeV). Thus the measurement of BR(h -> 



= 0(100 GeV) and 
t + t~) can provide 



a distinction between the SM and the MSSM even for relatively large M A and large nig if 
also and tan f3 are sufficiently large. The branching ratio BR(h — > bb) on the other hand, 
changes only by a few per cent for these parameters, so that this change would be much 
harder to measure. 
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Figure 7: BR(h — > t + t ) is shown as a function of m g . The left plot corresponds to three 
different values o/tan/3 and /i = ±100 GeV. The right plot corresponds to tan/3 = 40 and 
the two scenarios M A = 100 GeV,^ = ±100 GeV and M A = 300 GeV,/x = -400 GeV. 
The other parameters in both plots are nig = 500 GeV, X t = 500 GeV, M A = 100 GeV, 
M 2 = 500 GeV, At, = A t . The Higgs boson masses given in the legend are averaged masses 
over the interval < m g < 1000 GeV where the variation of Mh is about ±1 GeV for 
the above considered parameters. The horizontal lines represent the SM result using the 
corresponding Higgs boson masses. 



3.3 The a e ff-approximation 

In this section we investigate the quality of the a c fr-approximation. In Fig. |§| we display 
the relative difference between the full result ([34] ) and the a e fj result, where the external 
momentum of the Higgs self-energies has been neglected, see eq. (|35|) . The relative difference 
AT(h -> bb) = (T lull (h -> bb) - T a * s (h -> bb))/T lull (h -> bb) is shown as a function of M A 
for niq = 1000 GeV and for three values of tan/3 in the no mixing and the maximal mixing 
scenario. Large deviations occur only in the region 100 GeV < M A < 150 GeV, especially 
for large tan/3. In this region of parameter space the values of M^ and Mh are very close 
to each other. This results in a high sensitivity to small deviations in the Higgs boson 
self-energies entering the Higgs-boson mass matrix (^), (|T~5|) . Another source of differences 
between the full and the approximate calculation is the threshold M A = 2 m t = 350 GeV in 
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Figure 8: AT(h -> bb) = (T iuU (h -> bb) - T a ^(h -> bb))/T iull (h -> bb) is shown as a function 
of Ma for three values of tan (3. The QED, gluon- and gluino- contributions are neglected 
here. The other parameters are fi = —100 GeV, M 2 = nig, rrig = 500 GeV, Ab = A t , 
tan/3 = 3,20,40. The results are given in the minimal- and the maximal-mixing scenario. 
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Figure 9: The relative difference (sina c fr — sina c g(approx))/ sina e fl (see eqs. ffiBJ and ftSBj)) 
is shown as a function of Ma for three values of tan f3 in the no mixing and the maximal 
mixing scenario. The other parameters are chosen as in Fig. 



the one-loop contribution, originating from the top-loop diagram in the A self-energy and in 
the AZ self-energy (see Ref. J|). Here the deviation can amount up to 6%. 

In Fig. |9] we compare the a e s result ( |35D with the a e fr(approx) result (0), where the 
Higgs boson self-energies have been approximated by the compact analytical expression 
obtained in Ref. p2|. Fig. displays the relative difference in the effective mixing angles, 



18 



(sina eff — sina e fj(approx))/ sina eff . Via eq. fl26|) sina e fj directly determines the decay width 
T(h — > bb). The result is shown for rrig = 1000 GeV, for three values of tan/3 in the minimal 
and the maximal mixing scenario. Apart from the region around Ma ~ 120 GeV (compare 
Fig. D) both effective angles agree better than 3% with each other. 

Concerning the comparison of the a e fj-approximations in terms of Mh (which is not 
plotted here), due to the neglected external momentum or the neglected subdominant one- 
or two-loop terms, Mh receives a slight shift. Besides this kinematical effect, the decay rate 
is approximated rather well for most of the Mh values: independently of rrig, the differences 
stay mostly below 2-4%, for the no-mixing case as well as for the maximal-mixing case. Only 
at the endpoints of the spectrum, due to the different Higgs-boson mass values, the difference 
is not negligible. 



4 Comparison with the renormalization group approach 

In order to compare our results with those obtained within the renormalization-group- 
improved effective field theory approach |5],|| (in the following, for brevity reasons, denoted 
as RG approach), we made use of the program HDECAY In Fig. [II] we show the decay 
rate T(h — > bb) as a function of Mh (left part) and as a function of Ma (right part). We 
compare the RG result of HDECAY, where the gluino-exchange contribution is not yet im- 
plemented, with the FD result without and with gluino correction. In order to trace back 
the source of deviations we also show in the lower part of Fig. |10| the sine of the effective 
mixing angle a e s which enters the decay rate T(h — > bb) quadratically (see eqs. (^6|) and 
([ill)). In order to evaluate a e s we have neglected the external momentum (see eq. (|3~5"D ). As 
a typical soft SUSY-breaking scale we have chosen rrig = 500 GeV in Fig. |10[ Since the FD 
and the RG result have been obtained in different renormalization schemes, the entries of 



the i-mass matrix have a different meaning starting from two-loop order [^,|33],|34]]. For the 
no-mixing case we have set X t = in both approaches, whereas the maximal-mixing case is 
defined via X t = 2 rrig for the FD approach and X t = \/6rrig in the RG approach p2 |. 



In the comparison of the decay rates (and restricting ourselves to the results with ne- 
glected gluino exchange), for a given Mh we find deviations for low tan/3 up to O(10%) 
and agreement better than 4% in the large tan j3 scenario. The main part of the deviations 
can be attributed to the deviations in a e g which modify T(h — > bb) (see the lower part of 
Fig. |II]). As a general feature, larger deviations arise at the endpoints of the Mh spectrum 
due to the fact that for the same value of Ma different values for Mh are obtained in the 
FD and the RG approach, as shown in the left part of Fig. [TO. The plots for the tan f3 = 40 



scenario show again a sizable effect of the gluino correction. 

In Figs. and Fig. |l| the relative differences AT(h -> bb) = (T FD (h -> bb) - T RG (h -> 
bb))/T FD (h — > bb) and (sin a™— sin a^p) / sin a™ are shown as a function of Ma- Comparing 
a e ff as a function of Ma the agreement is relatively good; sizable deviations larger than 10% 



are found only in the regions where the curves in the right part of Fig. |10| have a steep 
slope. This can give rise to large deviations up to 50% in T(h — > bb) in terms of Ma (for 
100 GeV < Ma ^ 150 GeV), as displayed in the right parts of Figs. |TTJ and ^ (large tan f3). 
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Figure 10: F(h — > bb) is shown as a function of M h in the left part and as a function of Ma 
in the right part. The results of the RG approach (without gluino contribution) are compared 
with the Feynman diagrammatic results (without and with gluino contribution.) The other 
parameters are /i = —100 GeV, M 2 = rriq, rrig = 500 GeV, A b = A t . In each plot the 
no-mixing and the maximal-mixing scenarios are shown. In the lower part the sine of the 
effective mixing angle is shown. 
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Figure 11: AT{h -> bb) = (T FD {h -> 66)-r RG (/i -> 66))/r FD (/i -> 66) zs sftowi as a function 
of Ma for no mixing and maximal mixing. The gluino- contributions are neglected here. The 
other parameters are \i = —100 GeV, M 2 = m 9 -, m g = 500 GeV, Af, — A t . The results are 
given in the low tan/3 and in the high tan/3 scenario. 
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Figure 12: The relative difference (sin a™ — sina^)/ sin a™ is shown as a function of 
Ma for two values of tan f3 in the no mixing and the maximal mixing scenario. The other 
parameters are chosen as in Fig. [7J. 
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5 Conclusions 



Using the Feynman-diagrammatic approach for the Higgs-boson propagator corrections, in- 
cluding besides the full one-loop result also the dominant two-loop corrections, we have cal- 
culated the decay rates and branching ratios for the decays h — > bb, h — > cc and h — > t + t~. 
We have included the one-loop QED and QCD corrections, where the latter are due to gluon 
and gluino-exchange contributions. The gluino-exchange corrections have been neglected in 
most of the previous phenomenological analyses. In the present analysis, only the purely 
weak 0(a) (process specific) vertex corrections, shown to contribute less than 1% for most 
parts of the MSSM parameter space, have been neglected. 

We have shown analytically that the full set of Higgs-boson propagator corrections for 
vanishing external momentum can be absorbed into the effective mixing angle, a e g, in the 
neutral CP-even Higgs sector, appearing in the Higgs-boson fermion couplings. 

Numerically we have shown that, compared in terms of M h , the two-loop contributions 
to the Higgs-boson propagator corrections lead to a sizable decrease for T(h — > bb) and 
F(h — > t + t~), whereas F(h — > cc) can be increased, although it stays relatively small 
for all sets of parameters we have investigated. A sizable difference in all analyses from 
the one- to the two-loop calculation arises from the kinematical effect that at the two- 
loop level lower values for M h are obtained compared to the one-loop case, thus leading to 
deviations at the endpoints of the shown M h spectra. For most parts of the MSSM parameter 
space the a e fr-approximation reproduces the full result better than 3%. The gluino-exchange 
contribution to T(h — > bb) has been shown to be sizable for large tan/5. This correction 
increases with rising rrig and decreases with rising rriq and Ma- It is positive (negative) for 
negative (positive) /i. In the tan/5 = 40 scenario for small rriq it nearly compensate the 
gluon-exchange correction. 

The effect of the Higgs-propagator corrections on the branching ratios BR(h — > //) is 
relatively small for small values of Mh, while sizable effects are possible in the experimentally 
favored region of Mh- The branching ratio BR(h — > t + t~) can receive large corrections up to 
50% due to SUSY-QCD corrections to F(h — > bb). In some parameter regions the effective hbb 
coupling can become very small, and BR(h — > bb) can approach zero. In these parameter 
regions the two-loop contributions as well as the effect of the momentum dependence of 
the one-loop contributions are particularly important. It should be noted that in order to 
determine BR(h — > bb) very precisely in these paremeter regions, also a resummation of 
the leading terms as well as the inclusion of the complete 0(a) vertex corrections will be 
necessary. 

We have compared our diagrammatic result for Y(h — > bb) with the result obtained 
within the renormalization-group-improved effective field theory approach. For the low tan (3 
scenario, compared in terms of Mh, we find deviations up to 0(10%). In the large tan/5 
scenario the agreement is better than 4%. Compared in terms of Ma, however, differences 
of up to 50% are possible in the region 100 GeV < Ma ^ 150 GeV. The main part of the 
deviations can be attributed to the differences in the effective mixing angle a e ff- 
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